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dij 要 为 了 从 微观 分 子 相互 作用 层面 研究 三 维 打印 骨 支 架 工 艺 中 的 粉末 烙 结 机 理 及 本 质 , 本 文采 用 分 子 动力 学 的 模拟 仿真 
方法 , 分 别 从 内 聚 能 密度 结合 能 、 对 关联 函数 、 力 学 性 能 等 方面 对 目前 应 用 较 多 的 PVP、PAM、PVA 三 种 粘 结 剂 的 性 能 进行 了 
研究 , 并 将 所 得 结果 进行 了 分 析 和 比较 。 仿 真 结果 表明 , 三 种 粘 结 剂 与 HA 相互 作用 模型 的 界面 结合 能 的 大 小 关系 与 粘 结 剂 
本 身 的 内 聚 能 密度 大 小 关系 一 致 , 即 PAM > PVA > PVP; 粘 结 剂 高 分 子 与 羟基 础 灰 石 (Hydroxyapatite, HA) 的 对 关联 函数 分 
析 表 明 , 粘 结 剂 与 HA 发 生 相 互 作 用 主要 是 通过 高 聚 物 中 的 极 性 官能 基 团 与 HA 中 的 Ca 原子、 羟基 发 生 作用 形成 离子 键 、 氢 
键 , 且 离 子 键 作 用 强度 较 大 ; 此 外 , 三 种 相互 作用 模型 各 个 方向 的 力学 性 能 较 单一 HA 有 所 降低 , 且 相 互 作用 模型 力学 性 能 
的 优 务 关系 为 PVA/HA > PAM/HA > PVP/HA, 这 一 结论 与 结合 能 的 大 小 关系 并 不 完全 一 致 , 这 说 明 相 互 作 用 模型 的 力学 性 
能 与 粘 结 剂 的 粘性 并 不 存在 特定 的 内 在 关系 。 
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ABSTRACT In order to understand the bonding mechanism of hydroxyapatite (HA) particles for the 3D 
printedbone scaffolds with binders, the performance of three commercial binders i.e.PVP, PAM andPVA- 
was studied by means of molecular dynamics simulationin terms of cohesive energy density, binding en- 
ergy and pair correlation function g(r), as well as mechanical properties. The results revealed that the re- 
lationship of the binding energies between the HA surface with the three binders is consistent with their 
cohesive energy densities, i.e. PAM » PVA » PVP. The analysis of g(r) indicated that the interfacial inter- 
actionof HA and binders could mainly be attributed to the ionic bonds and hydrogen bonds which formed 
between the polar atoms, functional groups in binder polymer and the Ca, -OH in HA, and the strength of 
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ionic bonds is larger. TheYoung's modulus for the three interaction types of binders/HA can be ranked as 
the following sequence: PVA/HA > PAM/HA »PVP/HA, which are all inferior to that of the single HA. This 
conclusion is not completely consistent with the ranking of the relevant binding energies, which means 
that there is no specific intrinsic relation between the mechanical properties of the three binding types 


and the related viscidity of binders. 


KEY WORDS inorganic non-metallic materials, hydroxyapatite, molecular dynamics, bonding mecha- 


nism, mechanical properities 


羟基 础 灰 石 是 动物 人体 牙 齿 和 骨骼 的 主要 无 
TL ik pk 41", 该 材料 的 结构 类 似 于 天 然 骨 组 织 , 有 着 
良好 的 机 械 性 能 、 生 物 活性 、 生 物 相 容 性 、 骨 传导 性 ， 
无 毒性 、 化 学 性 质 稳 定 , 植 入 体内 不 易 引 起 炎症 和 免 
ERMA, 可 以 作为 骨 移 植 材料 的 基本 组 成 成 分 , 是 
一 种 最 具 潜 力 的 人 体 骨 组 织 蔡 换 修复 材料 , 在 骨 组 
织 工程 研究 中 得 到 了 广泛 的 应 用 。 

目前 羟基 磷 灰 石 骨 支架 材料 的 制备 还 都 是 采用 
较为 传统 的 方法 , 不 仅 操 作 复 杂 、 效 率 低 , 还 无 法 满 
足 病 患 对 骨 支 架 材 料 的 个 性 化 需求 。 随 着 先进 制造 
技术 的 不 断 发 展 , 出 现 了 基于 计算 机 辅助 设计 和 制 
造 的 粉末 粘 结 三 维 打印 快速 成 型 技术 (Three-Dimen- 
sional printing, 3DP), 该 技术 通过 基于 微 滴 喷射 原理 
的 三 维 打印 设备 喷射 液 滴 沉 积 于 粉 床 而 使 粉 材 粘 结 
成 型 具有 操作 简单 、 可 靠 性 高 、 成 型 速度 快 以 及 
环境 适应 性 好 等 优点 , 在 产品 的 开发 及 小 批量 生产 
中 具有 广阔 的 应 用 前 景 "特别 是 在 骨 文 架 制 备 
方面 已 得 到 了 成 熟 的 应 用 1。 而 在 粉末 粘 结 三 维 
打印 骨 支 架 的 制备 过 程 中 , 液 滴 与 粉末 的 相互 作用 
决定 了 骨 支 架 制 件 的 微观 结构 和 宏观 机 械 性 能 ， 
因此 在 HA 肯 支 架 的 制备 过 程 中 烙 结 剂 的 合理 选取 
显得 尤为 关键 。 粘 结 剂 不 仅 需 要 满足 对 材料 微观 结 
构 、 机 械 性 能 要 求 , 还 必须 是 安全 无 毒 、 对 HA 无 污 
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cohol, PVA) 三 种 胶水 , 这 3 种 胶水 具有 优良 的 生物 
惰性 、 生 物 相 容 性 和 良好 的 综合 性 能 , 能 很 好 地 满足 
粘 结 剂 对 HA 无 污染 的 要 求 。 

PAM、PVA、PVP 三 种 高 分 子 材 料 在 很 多 研究 中 
都 是 作为 水 凝 胶 、 薄 膜 主 体 材 料 来 进行 研究 的 , 对 其 
粘 结 行为 的 研究 鲜 有 报道 。 所 以 本 文 欲 通过 分 子 动 
力学 的 基本 理论 方法 对 PVP、PAM 以 及 PVA 在 HA 
界面 的 粘 结 行为 进行 深入 研究 , 从 分 子 相 互 作用 的 
角度 阐述 胶水 在 HA 表面 的 粘 结 机 理 及 本 质 , 并 将 
材料 的 各 项 性 能 参数 进行 相互 比较 , 分 析 、 总 结 3 种 
胶水 的 部 分 性 能 , 为 快速 成 型 骨 支架 制备 工艺 的 深 
入 研究 以 及 胶水 的 选取 提供 可 靠 的 理论 依据 。 

1 模型 的 建立 与 模拟 方法 

1.1 模型 的 构建 

首先 , 根据 PVP、PAM、PVA 的 分 子 结构 式 , 通过 
Materials Studio 软件 包 (Accelrys, San Diego, CA) 中 
的 Visualizer 模块 分 别 构建 PVP、PAM、PVA 高 分 子 
链 , 端 基 碳 原子 加 氧 以 达到 饱和 。 为 了 使 所 建 模型 
之 间 形 成 对 照 , PVP、PAM、PVA 高 分 子 链 聚 合 度 的 
选取 应 使 得 其 分 子 量 接近 相等 , 同时 三 种 粘 结 剂 高 
分 子 聚 合 度 ( 分 子 链 中 所 含 重 复 单 体 数目 ) 的 选取 必 
须 满足 高 分 子 聚合 度 的 最 小 要 求 , 即 能 够 代表 高 分 


m 
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染 , 目前 应 用 较 多 的 是 医用 生物 胶水 (a - 氰 基 丙 稀 
酸 酯 类 ), 但 医用 胶水 的 成 本 过 高 、 容 易 破坏 肯 支 架 
的 部 分 微观 结构 , 并 不 适合 在 骨 文 架 材 料 生产 中 大 
使 用 。 所 以 普遍 采用 的 方法 是 先 用 其 它 胶水 来 制 


子 的 最 小 重复 单元 数 。 本 文通 过 计算 , 在 满足 上 述 
条 件 下 , 将 PVP、PAM、PVA 三 种 高 聚 物 的 聚合 度 分 
HJA 20,31,50. XH Smart Minimizer 方 法 对 其 
进行 优化 , 然后 利用 Amorphous Cell 模块 分 别 构建 


gm jl 


rx RBS ABE, 然后 用 医用 胶水 对 所 制备 的 骨 文 
架 和 雏形 进 行 加固 处 理 , 以 提高 其 机 械 性 能 。 目 前 使 
用 最 多 的 是 聚 乙烯 毗 咯 煤 酮 Polyvinylpyrrolidone， 
PVP)、 聚 丙烯 酰胺 OPAM) 以 及 聚 乙烯 醇 (polyvinyl al- 
表 1 REWER 


Table 1 Detail parameters of 


三 种 高 分 子 材料 的 无 定型 唱 胞 模型 。 为 了 减少 尺寸 
效应 , 而 又 不 至 于 使 计算 量 过 大 , 每 个 无 定型 唱 胞 含 2 
条 高 聚 物 链 , 所 对 应 的 唱 胞 参数 如 表 1 所 示 ; 此 外 , f 
据 HA 的 品格 参数 , 构建 HA 的 单元 体 唱 胞 模型 , 所 建 


型 晶 胞 的 详细 参数 


amorphous cells for polymers 


Bite Number of Number of Number of Initial density/g* Finial density 
repeat units chains atoms cm? /gXcm? 
PVP 20 2 684 0.6 1.201 
PAM 31 2 624 0.6 1.316 
PVA 50 2 704 0.6 1.265 
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单元 体 晶 胞 模型 的 晶 胞 参数 为 : a— 50.9432 nm, c= 
0.6881 nm, a=B=90", y=120°"". F HA 是 一 种 各 
各 异性 材料 , 本 文 只 选取 其 原子 密度 最 大 的 一 个 平 
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党 Od 30 着 


了 三 种 高 分 子 材料 与 HA(110) 的 界面 相互 作用 的 平 
衡 结 构 模 型 , 如 图 1 所 示 。 其 中 紫色 为 高 聚 物 的 碳 
主 链 , 灰色 为 体系 中 主 链 碳 原子 外 的 碳 原 子 , 红色 为 


面 构建 超级 晶 胞 模型 9, 即 最 具 代 表 性 的 HA(110) 平 
mo KE, 截取 HA 单 唱 胞 (110) 平 面 , 重新 构建 超级 晶 
胞 , 所 建 超级 晶 胞 的 晶 胞 参数 为 a=3.2 nm, 572.0 nm, 
c=1.1 nm, a=B=y=90"。 最 后 通过 Build Layers 模块 
构建 高 聚 物 与 HA(110) 的 相互 作用 模型 。 
1.2 MD 模拟 细节 

无 论 是 粘 结 剂 的 无 定型 晶 胞 , 还 是 其 与 HA 界 


原子 氧 , 乳白 色 为 氧 原 子 , 蓝 色 为 所 原子 , E GS E 
原子 , 绿色 为 钙 原子 。 

2 结果 与 讨论 
2.1 内 聚 能 密度 与 分 子 链 聚合 度 的 选取 
真实 的 高 分 子 链 所 含 原 子 个 数 成 干 上 万 , 在 模 
拟 过程 中 构建 一 条 真实 的 高 分 子 链 是 不 切实 际 的 ， 


面 的 相互 作用 模型 , 在 进行 数值 分 析 之 前 , 都 必须 
使 其 达到 平衡 状态 , 即 能 量 最 小 状态 , 这 关系 到 模 
拟 结果 的 准确 性 。 分 子 动力 学 模拟 有 具体 包括 以 下 2 
个 段 : 

(1) 平衡 模拟 阶段 , 该 过 程 首 先 选 用 的 是 NVT 
REK, 温度 设 为 298 K, 模拟 时 间 设 置 为 60 ps, 以 使 
体系 达到 松弛 。 而 后 再 在 NPT 系 综 下 , 将 压强 设置 
23 1 bar, 温度 设置 为 298 K, 模拟 时 间 为 100 ps, 以 使 
体系 达到 最 终 的 平衡 构象 。 

Q) 分 析 模 拟 阶段 : 在 平衡 模拟 的 基础 上 , 通过 
NVT 系 综 分 子 动 力学 模拟 , 模拟 时 间 设 置 为 80 ps, 
并 将 其 中 的 后 30 ps 的 轨迹 文件 用 于 数值 分 析 。 

分 子 动力 学 过 程 都 是 在 COMPASS 力 场 下 进行 
的 , 3€ H] Andeson TH 2 38 ^, 用 Velocity Verlet 积分 法 


这 是 由 计算 能 力 的 客观 因素 所 决定 的 。 因 此 , 高 分 
子 链 聚 合 度 的 选取 是 模拟 成 功 的 关键 , 高 分 子 链 所 
含 单元 体 个 数 既 不 能 超出 计算 量 的 上 限 , 还 必须 能 
够 代表 高 分 子 的 性 质 。 内 聚 能 密度 (Cohesive Ener- 
gy Density, CED) 是 表征 高 分 子 材料 的 一 个 重要 参 
数 , 是 对 物质 分 子 间 相互 作用 力 大 小 的 度量 , 几 与 物 
质 间 相互 作用 有 关 的 性 质 和 物性 属性 与 内 聚 能 密度 
都 具有 内 在 联系 , 如 物质 的 溶解 性 、 相 容 性 、 汽 化 热 
等 , 其 表达 式 如 下 : 
E 


CED = mE (1) 


AP Ew 是 指 内 聚 能 , v 为 物质 的 量 混合 体积 。 而 重 
复 单元 体 个 数 与 高 聚 物 的 内 聚 能 存在 一 定 的 关系 ， 
即 随 着 重复 单元 体 个 数 的 增加 , 高 分 子 链 的 CED 趋 


求解 牛顿 运动 方程 , 原子 初始 速度 按 Boltzmann 随 
机 分 布 方法 确定 。 计 算 采 用 的 时 间 步 长 为 1 fs, 模 
拟 过 程 中 每 500 步 记录 一 次 体系 的 轨迹 。 模 拟 过 程 
中 Vander Waals fll Coulomb 分 别 按 Atom based" 和 
Ewald" 7; 15 ib $6, 截断 半径 为 0.95 nm, 样 条 宽度 取 
0.1 nm, 缓冲 宽度 取 0.05 nm, 温度 设置 为 298 K. 
1.3 平衡 结构 
体系 的 平衡 可 由 温度 和 能 量 同时 到 达 平 衡 予 以 
确定 。 通 常 当 温度 和 能 量 在 5%~10% 范 围 内 波动 即 
可 认为 体系 已 达到 平衡 。 本 研究 的 所 有 计算 都 是 在 
体系 达到 平衡 的 条 件 下 进行 的 。 经 MD 模拟 , 得 到 
(a) 


FE, 所 以 可 以 根据 高 分 子 链 聚 合 度 与 内 聚 能 密 
度 的 关系 来 确定 代表 真实 高 分 子 链 的 最 小 链 长 。 
为 了 选取 合适 的 高 分 子 聚合 度 , 本 文 构建 了 含 
不 同 聚 合 度 的 PVP、PAM、PVA 分 子 链 模 型 , 并 对 其 
无 定型 晶 胞 进行 了 MD 模拟 , 计算 得 到 了 不 同 聚 合 
度 高 分 子 链 的 内 聚 能 密度 , 图 2 为 内 聚 能 密度 与 高 
分 子 聚 合 度 的 关系 。 
通过 内 聚 能 密度 与 聚合 度 的 关系 图 , 可 见 PVP、 
PAM 的 聚合 度 超过 20 后 , 内 聚 能 密度 就 趋 于 稳定 ， 
而 PVA 则 要 超过 30 才能 趋 于 稳定 。 而 本 文 为 了 对 
三 种 材料 进行 对 比 , 同时 权衡 计算 能 力 , 最 终 将 


图 1 三 种 粘 结 剂 与 HA 718 


相互 作用 模型 的 平衡 结构 


Fig.1 Equilibrium structures of binder/HA interfacial interaction models. (a) PVP/HA(110) interaction 
model; (b) PAM/HA(110) interaction model; (c) PVA/HA(110) interaction model 
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为 了 进一步 前 述 粘 结 剂 与 HA 界面 间 粘 结 机 理 
及 本 质 , 本 文 对 三 种 粘 结 剂 与 HA 界面 相互 作用 体 
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S o8 t Z RAM 1 4& PVP f A 20 MIKE (-C-O), 1 条 PVA 链 含 有 
z 07H 50 个 羟基 (-OH), T 1 4& PAM 链 则 含有 31 个 酰胺 基 
B oal CNH2).31 433&(-C-O). 
ZI 2.3 对 关联 函数 分 析 
L 
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10 20 30 40 50 60 
The number of repeat units 

2 内 聚 能 密度 随 聚 合 度 的 变化 

Fig.2 CED vs repeat unit number 


o 


o 


PVP, PAM, PVA 的 聚合 度 取 为 20、31、50, 即 单 链 的 
分 子 量 相 近 。 此 外 通过 比较 三 种 高 聚 物 稳定 时 的 
CED, 其 大 小 关系 为 PAM>PVA>PVP; mi CED 的 大 
小 可 以 反映 物质 的 稳定 性 以 及 其 内 部 分 子 间 相 互 作 
用 的 强 弱 , 即 内 聚 能 密度 越 大 物质 的 稳定 性 、 内 部 相 
互 作用 力也 就 越 大 , 所 以 高 聚 物体 系 的 稳定 性 及 其 
内 部 相互 作用 的 强 弱 关 系 为 PAM>PVA> PVP。 


2.2 结合 能 


高 聚 物 与 HA(110) 界 面 的 结合 能 是 高 聚 物 分 子 

与 HA 界面 分 子 相互 作用 力 大 小 的 一 个 量化 , 可 以 

反映 出 高 聚 物 粘 结 剂 与 HA 表面 的 粘 结 作 用 力 大 
小 , 是 粘 结 剂 的 一 个 重要 参数 。 其 表达 式 如 下 : 

Epira Eum rc E Pula E ua) (2) 

式 中 , Esa Erme ^ Enno 分 别 表 示 高 聚 物 与 HA 


相互 作用 体系 、 高 聚 物 .HA(110) 的 能 量 。 通 过 对 
PVP/HA(110)、PAM/HA(110)、PVA/HA(110) 的 相互 


系 的 平衡 结构 进行 对 关联 函数 (Pair Correlation 
Function, PCF) 分 析 。 对 关联 函数 ge(D 可 表示 在 距离 
某 一 指定 参考 原子 A 为 + 处 , B 原子 出 现 的 概率 , 其 
表达 式 具 体 如 下 。 


ga-slr)=( 


ng Ns 

iu y (C77) (3) 
式 中 ne 表示 在 距离 A 原子 半径 为 r 的 范围 内 B 原子 
的 数目 , Ns 表示 整个 体系 中 B 原 子 的 总 数 , v 表示 整 
个 体系 的 体积 。 

在 进行 对 关联 函数 的 分 析 之 前 , 首先 必须 对 相 
互 作 用 体系 中 的 各 类 原子 进行 标记 , PVP 分 子 链 中 
的 原子 类 型 有 主 链 C 原子 .与 主 链 C 原子 直接 相连 
WHET, N 原子 .与 N 原子 构成 五 元 环 的 C 原子 以 
J& FE (-C—0)O 原子 , 其 中 只 需 对 可 能 与 周围 原子 
发 生 较 强 相 互 作用 的 原子 进行 标记 , 所 以 将 与 主 链 
C 原子 直接 相连 的 互 原 子 `N JEU EC C-0)0 原 
子 依次 标记 为 HPVP) N(PVP), O(- C20); 类 似 的 ， 
在 PAM 中 将 与 主 链 C 原子 直接 相连 的 五 原子 、 装 基 
(-C=0)O 原子 .酰胺 基 (-NHJN 原子 、 栈 胺 基 瑞 原子 
标记 为 H(PAM), O(- C20). N(- NH), H(- NH); 在 
PVA 分 子 链 中 , 将 羟基 (-OH)JO Fi T  X83E H CT. 5 
主 链 C 直 接 相 连 的 互 原子 依次 标记 为 O(C-OH) 、HC- 


作用 平衡 体系 进行 MD 模拟 计算 , 所 得 到 的 结合 能 
如 表 2 Bra: 
由 表 2 数据 可 见 , PVP, PAM, PVA 三 种 粘 结 剂 
高 分 子 材料 在 HA 表面 的 结合 能 大 小 关系 为 PAM> 
PVA>PVP。 说 明 这 三 种 粘 结 剂 与 HA 界面 粘 合力 的 
大 小 关系 为 PAM>PVA>PVP, 即 PAM 对 HA 界面 的 
粘 结 力 相 对 最 强 。 而 组 分 间 结 合 能 的 大 小 是 界面 分 
子 相 互 作用 力 的 反映 , 结合 能 越 大 , 产生 的 分 子 间 相 
互 作用 力 越 大 , 这 与 内 聚 能 密度 的 大 小 关系 相 一 
致 。 导 致 上 述 结果 的 最 根本 原因 是 因为 相同 质量 的 


表 2 高 聚 物 与 HA(110) 界 本 


OH), H(PVA); 将 HA 中 的 Ca 原子、 磷酸 根 (-PO7) 中 
的 O 原子 `H 原 子 和 羟基 (COH) 中 的 O 原子 , 依次 表 
7R 73 Ca. O(HA). H(HA), OHA); 本 文 分 别 对 上 述 
三 种 高 聚 物 所 标记 的 原子 与 HA 中 所 标记 的 原子 
之 间 进 行 了 对 关联 函数 分 析 ( 图 3-4), 此 外 还 对 三 种 
粘 结 剂 高 分 子 C 主 链 在 HA 界面 的 分 布 进 行 了 分 析 
(图 3)。 图 4 只 给 出 了 存在 重要 相互 作用 的 gaor 
图 3 为 三 种 粘 结 剂 高 分 子 C 主 链 与 HA 界面 的 
对 关联 函数 , h d n] JU PVP. PAM. PVA 的 C 主 链 所 


i 相互 作用 体系 的 结合 能 


Table 2 The binding energies between polymers and HA(110) surface (kJ/mol) 


System Fa Fuse: Eus Eina 
PVP/HA 398067331 36773 39917761 147803 
PAM/HA 39706547 58660 39917652 269765 
PVA/HA 39744011 49373 39917778 223140 
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(a) 
— OCCc=O)-C 
14 — C(PVP)-HA -1 = eed - 
12 — CPAM HA — N(PVP)-Ca 
SP RAUS sl 一 NIPVP)-HHA) 
gis " — H(PVP)-O(HA2) 
S og S 6 
0.6 4 
0.4 
2 - 
0.2 
0 0 
0 02 0.4 0.6 0.8 0 
ri nm 10 
(b) a 
图 3 主 链 C 原子 与 HA 界面 的 对 关联 函数 a — o ESHE 
Fig.3 The PCFs between the carbon backbone and HA in- — N(-NH2)-Ca 
terface — N(-NH2)-H(HA) 
c BT —  H(NH2)-O(HA2) 
D 


对 应 的 对 关联 函数 值 的 大 小 关系 为 PAM>PVA> 
PVP, rfj C 主 链 可 以 在 一 定 程 度 上 代表 整个 高 分 子 ， 
所 以 函数 值 的 大 小 可 以 反映 出 粘 结 剂 与 HA 界面 粘 
合力 的 大 小 , 即 相 互 作用 力 越 大 , 其 在 距离 HA E I 
同一 距离 7 处 的 函数 值 也 越 大 。 这 一 结论 与 粘 结 剂 
与 HA 相互 作用 结合 能 的 大 小 关系 是 一 致 的 , 进 一 
步 验 证 了 计算 结果 的 正确 性 。 

分 子 间 相 互 作用 力 一 般 分 为 化 学 键 、 氧 键 作 用 
和 范 德 华 作用 , 对 关联 函数 图 中 0.35 nm 以 内 的 峰值 
主要 由 化 学 键 、 氧 键 构成 , 0.35 nm 以 外 的 则 主要 是 
库伦 、 范 德 华 相 互 作用 成 分 时。 图 4a 分 析 了 PVP/ 
HA 界面 相互 作用 模型 中 存在 重要 相互 作用 的 原子 
间 对 关联 函数 , PVP fh RP EEECC-O)O 原子 、 
N 原子 与 HA 中 Ca 原子 的 对 关联 函数 的 峰值 位 于 天 
0.22 nm 附近 , 比 Ca-O 离子 键 键 长 0.239 nm、Ca-N 离 
子 键 键 长 0.247 nm 均 小 , X HH PVP P S Exc O 原 
子 、N 原子 分 别 与 HA 中 的 Ca 原子 之 间 形 成 了 离子 
; 此 外 , PVP rP 3X O 原子 、N 原子 与 HIHA) 原 
对 关联 函数 的 峰值 位 于 770.21 nm 附近 , 说 明 它 们 


AA 


一 O(OH)-Ca 
8 上 —  O(-OH)-H(HA) 
一 H(OH-O(HA2) 
一 H(PVA)-O(HA2) 
一 H(PVA)-O(HA1) 


g) 


rinm 


4 AAEH A 界面 相互 作用 的 对 关联 函数 

Fig.4 PCFS of binder/HA interface interaction: (a) PCFS 
of PVP/HA interface interaction; (b) PCFS of PAM/ 
HA interface interaction; (c) PCFS of PVA/HA in- 
terface interaction 


之 间 主 要 是 通过 氧 键 发 生 相 互 作用 ; 而 HPVP) 与 O 
(HA2) 对 关联 函数 的 峰值 不 是 特别 明显 , 但 是 也 形 
成 了 少 部 分 氧 键 , 大 部 分 是 范 德 华 作用 ; 由 上 述 原子 
间 对 关联 函数 的 峰值 大 小 , 可 以 发 现 PVP 与 HA 界 
面相 互 作用 主要 是 靠 PVP 中 的 极 性 原子 官能团 与 
HA 界面 原子 形成 离子 键 \, 氧 键 , 且 离 子 键 的 作用 强 
度 最 大 。 

B| 4b J PAM/HA 界面 相互 作用 的 对 关联 函数 ， 
由 图 可 见 , PAM HP B5 d O 原子 酰胺 基 KN 原子 与 
HA 中 Ca 原子 间 的 对 关联 函数 峰值 分 布 在 天 0.22 nm 
Bir, 说 明 它 们 之 间 形 成 了 离子 键 ; 此 外 , xd O 原 
子 、 栈 胺 基 N 原 子 与 HA 中 羟基 五 的 对 关联 函数 峰 
值 也 在 氧 键 的 作用 范围 之 内 , 说 明 它 们 之 间 形 成 了 


Ec 


^ f fF Hd; 而 酰 胶 基 瑞 原 子 与 HA rn 83€ O 原子 间 
的 相互 作用 不 是 很 强 , 主要 是 发 生 范 德 华 作 用 ; 从 上 
述 原 子 间 对 关联 函数 的 峰值 大 小 分 析 可 以 看 出 
PAM 与 HA 界面 相互 作用 主要 是 通过 默 基 O 原子 、 
酰胺 基 N 原子 与 HA 界面 原子 形成 离子 键 、 氧 键 , Ho 
离子 键 作 用 强度 最 大 , 其 次 是 氧 键 作 用 。 

图 4c 为 PVA/HA 界面 相互 作用 的 对 关联 函数 ， 
由 图 可 见 PVA 中 羟基 (-OHJO 原子 与 HA 中 的 Ca Ji 
子 、 羟 基 瑟 原子 间 的 对 关联 函数 峰值 分 别处 于 产 
0.21~0.22 nm 附近 , 说 明 其 分 别 形成 了 离子 键 、 氢 键 
作用 ; 而 PVA 中 羟基 卫 原子 与 HA "PESE O 原子 间 
的 对 关联 函数 的 峰值 位 于 二 0.31 nm 附近 , 这 是 氧 键 
的 作用 范围 , 说 明 它 们 之 间 是 通过 氧 键 发 生 相 互 作 


-D 


ao 


202303.10700v1 


chinaXiv 


8 期 


pA 


Hj; PVA rp 5; xx SE C ET ERIE H AT 5 HA 


rp O 原子 间 的 相互 作用 很 微弱 ， 
比较 对 关联 函数 峰值 的 大 小 可 以 发 现 PVA 与 HA 界 


主要 是 范 德 华 作用 ; 


面相 互 作 


也 主要 是 通过 PVA 中 的 羟基 与 HA JR 


原子 形成 离子 键 、 氨 键 , 且 离 子 键 作用 强度 最 大 , AG 


所 以 通过 对 PVP、PAM、PVA 三 种 粘 结 剂 高 分 子 
与 HA 的 对 关联 函数 图 像 分 析 , 可 以 发 现 , 粘 结 剂 与 


HA 发 生 相 互 作 月 
基 团 与 HA 中 的 Ca 原子 \ 凑 基 发 生 作 用 形成 离子 键 、 


主要 是 通过 高 聚 物 中 的 极 性 官能 


AIE, 且 离 子 键 作 用 强度 最 大 , 氧 键 次 之 。 


2.4 力学 性 能 分 析 


肯 文 架 的 制备 过 程 中 , 力学 性 能 是 


剂 的 一 个 重要 指标 , 本 文 应 


= 选 


用 粘 结 
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论 


1. 通过 计算 比较 三 
E, 得 出 了 三 


' 粘 结 剂 稳 定时 的 内 聚 能 密 


材料 的 稳定 性 及 其 内 部 相互 作用 的 


强 弱 关系 为 PAM>PVA>PVP。 


2. 结合 能 的 大 小 是 界面 分 子 相互 作用 力 大 小 的 


反映 , 结合 能 越 大 , 分 子 间 相 互 作 


j 力 就 越 大 , 所 以 


由 结合 能 的 大 小 关系 得 出 了 三 种 粘 结 剂 与 HA 界面 
粘 合力 的 大 小 关系 , 即 为 PAM>PVA>PVP, ix 5 N% 
能 密度 的 大 小 关系 相 一 致 。 


3. 通过 对 PVP、PAM、PVA 三 种 粘 结 剂 高 分 子 与 


HA 的 对 关联 函数 分 析 , 可 以 发 现 ， 


离子 键 作 有 


件 中 的 Forcite 模块 , 对 三 种 粘 结 剂 与 HA 相互 作 月 


H 


^F f ER X. 六 2 三 个 方向 (QZ 为 作用 面 的 垂直 方 


图 $ 所 示 。 


由 图 5 不 难 发 现 三 


粘 结 剂 与 HA 相互 作 


句 ) 的 杨 氏 模 量 进行 了 MD 模拟 计算 , 其 计算 结果 如 


用 模 


型 三 个 方向 的 杨 氏 模 量 并 不 一 致 , 存在 一 定 差 
3r, 这 一 结论 与 HA 的 各 向 异性 相符 , 即 粘 结 剂 的 


添加 3 
的 相互 作用 模型 


没有 改变 材料 的 各 向 异性 ; 添加 粘 结 剂 后 


,其 各 个 方向 的 力学 性 能 较 


、 
HR 


HA 力学 性 能 均 有 所 下 降 , 表现 最 为 明显 的 是 Z 
方向 的 力学 性 能 , 这 主要 与 建立 相互 作用 模型 


时 作用 国 
方向 添加 粘 结 剂 


的 选取 有 关 , 即 本 文 是 在 HA 超级 唱 胞 Z 


; 比较 三 个 相互 作用 模型 的 力学 


性 能 , 可 以 发 现 其 力学 性 能 存在 一 定 的 优 劣 关系 ， 


力学 性 能 最 好 的 是 PVA/HA, 其 次 是 PAM/HA, RIR 


的 是 PVP/HA, 即 PVA>PAM>PVP, 这 一 结论 与 结合 
能 的 大 小 关系 并 不 完全 一 致 , 这 说 明 相互 作用 模 
型 的 力学 性 能 与 粘 结 剂 的 粘性 并 不 存在 特定 的 内 
在 关系 。 
70 
" EA PAUNA EZ 
B EC 
S f z ES PVA/HA p 
$ sol 7 多 HA 
2 d NS p 
il gm N y 
:oN DN wN 
sI UN, AN oN 
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图 $ 相互 作 


Fig.5 Young's modulus histograms of the interaction models 


模型 的 杨 氏 模 量 柱状 图 


日 强度 最 大 , 氧 键 次 之 。 


粘 结 剂 与 HA 发 


生 相互 作用 主要 是 通过 高 聚 物 中 的 极 性 官能 
HA 中 的 Ca 原子、 羟基 发 生 作 用 形成 离子 键 、 所 键 ， 


基 团 与 


单 
相 
的 


4. 三 种 相互 作用 模型 各 个 方向 的 力学 性 能 较 
一 HA 有 所 降低 , 且 Z 方 向 尤为 明显 。 对 比 三 个 
互 作用 模型 的 力学 性 能 , 其 粘 结 模型 力学 性 能 
优 劣 关 系 为 PVA>PAM>PVP, 这 一 结论 与 结合 能 


的 大 小 关系 并 不 完全 一 致 , 这 说 明 相 互 作用 模型 
的 力学 性 能 与 粘 结 剂 的 粘性 并 不 存在 特定 的 内 在 
关系 。 
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